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We study both mean-field and full quantum dynamics of symmetry-breaking transitions (SBTs) 
in a coupled two-component Bose-Einstein condensate. By controlling s-wave scattering lengths 
and coupling strength, it is possible to stimulate SBTs between normal and spontaneously polarized 
ground states. In static transitions, the probability maxima of full quantum ground states correspond 
to the mean-field ground states. In dynamical transitions, due to the vanishing of excitation gaps, 
the mean-field dynamics shows universal scalings obeying Kibble-Zurek mechanism. Both mean- 
field and full quantum defect modes appear as damped oscillations, but they appear at different 
critical points and undergo different oscillation regimes. The anomalous breakdown of mean-field 
dynamics induced by SBTs depends on the approaching direction. 
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In a many-body quantum system, spontaneous symme- 
try breaking (SSB) occurs if its mean-field (MF) states do 
not possess symmetry of its original many-body Hamil- 
tonian. Since atomic Bose-Einstein condensates (BECs) 
have long collision times, they are excellent candidates 
for testing intrinsic mechanisms of the SSB []]]. In par- 
ticular, the simultaneous realization of superfluidity and 
magnetism (or spin polarization) in spinor BECs is asso- 
ciated with the SSB related to both external and internal 
degrees of freedom. For spin-1 BECs, the experimental 
observation of SSB 0] has triggered several theoretical 
studies 0, 0|. For spin-^ (two-component) BECs, the 
phase separation has been observed experimentally [H], 
and their spatial SSB Q and spontaneous spin polariza- 
tion have been predicted theoretically. 

It is well known that, in slow processes, the gapped 
excitations over zero-temperature ground states obey 
Landau-Zener mechanisms. However, gapless excita- 
tions appear in almost all systems with SSB. Because of 
the gapless excitations, the adiabaticitybreaks down j§| 
and the generated defect modes 0, 0, 11 follow Kibble- 
Zurek (KZ) mechanisms [IB]. The current studies of 
dynamical mechanisms for SBTs focus on lattice sys- 
tems @, S [ll|, spin-l BECs 0, 0] and other many- 
body systems by employing either MF 0, EH or full 
quantum (FQ) [3, 0, [T3| theories. However, (i) few 
works compare the MF and FQ dynamical mechanisms 
near a critical point to explore regimes of correspon- 
dence, and (ii) the dynamical mechanisms for SBTs in 
two-component BECs are still not clearly understood. 

In this Letter, we analyze both MF and FQ dynam- 
ical mechanisms for SBTs in a coupled two-component 
BEC. For simplicity, we only consider the SSB related 
to internal degrees of freedom, i.e., two hyperfine spin 
states. In static transitions, the MF ground states cor- 
respond to the probability maxima of the FQ ground 



states. If the coupling strength increases from zero to a 
sufficiently large quantity, the MF ground states trans- 
fer from being spontaneously polarized (self-trapped) to 
non-polarized (normal) . Correspondingly, the two lowest 
FQ eigenstates transfer from being quasi-degenerated to 
non-degenerated. In dynamical transitions, due to the 
disappearance of gaped Bogoliubov excitations, the MF 
dynamics obeys an universal KZ mechanism and the de- 
fect modes are induced by dynamical instability. Due 
to the non-identity of the FQ and MF critical points, 
the MF breakdown induced by SBTs depends on the ap- 
proaching direction. The dynamical mechanism of SBTs 
connects with the quantum adiabaticity, which provides 
various applications in atomic physics, condensed matter 
physics and non-equilibrium dynamics, and particularly 
in adiabatic quantum computation flU ]. 

We consider a gaseous BEC of bosonic atoms which 
populate two hyperfine levels coupled by external fields. 
Assuming that the coupling fields are spatial uniform and 
the atom-atom interactions do not change the internal 
states, the system obeys a second quantized Hamilto- 
nian [3], 

H(t) = H +H int + H c (t), 
H ° = E h+^i-^ + V + e^A^r, 

Hint = \ E U i>< I *+«*^(r)* fc (r)^ (r)d 3 r, 

Hc(t) = J (*+(f)^(r) + *+(r)* T (r))dV. 

Here, m is the single-atom mass, V — ^moj 2 r 2 is the 
trapping potential, fl(t) > is the coupling strength, Cj 
are the hyperfine energies, $j(r) are the field operators, 
and Ujk = Ukj = Anffajk/m are interaction strengths of 
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atoms in states \j) and \k) parameterized by s-wave scat- 
tering lengths djk- To analyze dynamics only associated 
with internal degrees of freedom, we apply a single-mode 
approximation [1J|: $j(r) = bj<f>(r), that is, assume that 
all atoms occupy the same single-particle external state 
4>{r). Thus the Hamiltonian is simplified as follows: 
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where E Qj = / (jf(r) (- 
G jk = U jk JW)\ 4 d 3 f 



V(r) 



(f)(r)d 3 r and 



This simplification can suc- 
cessfully describe the system of uniform polarization. 
The Hamiltonian is equivalent to a Bose-Josephson junc- 
tion [H, 0, Il7j . which can be realized by a BEC in a 
double- well trap In contrast to the double- well sys- 
tems, in which negative charging energies Ec oc a s may 
cause spatial collapse, the two-component systems sup- 
port negative Ec oc (a-ff + aj; — 2ajjJ with all positive 
a-jk- 

In MF theory, by introducing variational trial states 
|$) as coherent states [3], ipj = ($|S+|$) and ipj = 

obey the classical Hamiltonian: 
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Rewriting ipj = y/Nj exp(i9j) in terms of the particle 
numbers Nj and phases 6j, we denote I = (N^ — iVj.)/2, 
L = (N i + N r )/2, A = E ]- — E ol + L(G n ~ G u ) and 
Ec = Gii + ^11 — 2G||. It has been predicted that 
spontaneous spin polarization, a type of SSB related 
to self-trapping and bi-stability flif . [20| . can appear 
when Ec < 0. For symmetric systems (A = 0) 
of negative Ec, the SBT occurs if \HQ/Ec\ varies 
from \hVt/E c \ > L to \Ml/E c \ < L. Correspond- 
ingly, the MF ground state changes from (ipi , ipl) — 
(VL, y/L) exp (iOi ) to 
where L 



WL - l s , y/L + l s ) exp (iO^, 
Below, we will focus on 



±y/L 2 - h 2 n 2 /E' c . 

the SBTs in symmetric systems of Ec < 0. 

First, let us analyze the MF-FQ correspondence for 
static SBTs. It is obvious that the FQ ground state 
appears as an SU(2) coherent state if Ec — 0, and 
the ground and first excited states become degenerate 
if fi = and Ec < 0. For other arbitrary parameters, 
our numerical results show that the FQ ground states are 
always symmetric with respect to I = and the proba- 
bility distributions change from single-hump shapes to 
double-hump ones when Ml/Ec passes the critical point 
Ml/Ec = —L. However, the first excited states are al- 
ways anti-symmetric with respect to I = and the prob- 
ability distributions retain double-hump shapes. 



The probability distributions and degeneracy proper- 
ties of the low-energy FQ states are reminiscent of a single 
quantum particle confined within a potential that varies 
from single-well to double-well configuration. By com- 
paring the FQ and MF ground states, it is easy to find 
that the MF ground states correspond to the probability 
maxima of the FQ ones. Because the quasi-degeneracy 
between the two lowest FQ states requires almost iden- 
tical probability distributions for these two states, the 
MF bifurcation point is not identical to the degener- 
acy point for the two lowest FQ eigenstates. In Fig. 1, 
we demnstrate the MF-FQ correspondence for the static 
SBT in a symmetric system of total particle number 
N = N t + 7V T = 100 and E c < 0. The supercritical 
pitchfork bifurcation of the MF stationary states signals 
the static SBT. Similar bifurcations have also been found 
in rotating BECs [2l[. Due to the appearance or dis- 
appearance of unstable or stable stationary states, the 
bifurcation could cause dynamical instability. 




Figure 1: The MF-FQ correspondence for static SBTs in a 
coupled two-component condensate of N = 100 and nega- 
tive E c . The FQ states are defined as |#) = ^fL_ L C{l)\l) 
and the probability distributions |C(7)| 2 versus Q,/Ec for the 
ground (first-excited) states are shown in the above (bottom) 
panel. The solid and dashed lines are stable and unstable MF 
stationary states, respectively. 

Now, let us consider the dynamical SBTs. In a dynam- 
ical process, there exist two characteristic time scales 
One time scale is the reaction time, T r (t) — h/A g (t), 
which characterizes how fast the system follows eigen- 
states of its instantaneous Hamiltonian. Here, A 9 (t) 
is the instantaneous excitation gap over the ground 
state. The other timescale is the transition time, r t (t) — 
A g (t)/ \dA g (t)/dt\, which tells us how fast the system 
is driven. If and only if r r < r t , the system undergoes 
adiabatic evolution. 
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To obtain the excitation modes over MF ground states, 
we perform a Bogoliubov analysis. Apart from a trivial 
gapless mode, there exists a gapped mode: 



f h^/ttjtt 



tt? 



), for tt > tt c 

for tt < tt c , 



(3) 



where tt c = \EcL/h\. Obviously, the gap A g gradually 
vanishes when tt approaches the critical point tt c . 

To analyse the dynamical mechanism, we suppose 
tt(t) = tt c {\±t/T q ) = fi c ± f3t, where r q denotes 
the quenching time. The relative coupling, e = 
|(fi(i) — fi c ) /f2 c | = \t\/T q , corresponds to the relative 
temperature in KZ theory [IB]- Due to the gap A g — > 
when \tt — tt c \ — > 0, the adiabaticity breaks down when 
tt —y tt c and then revives after tt passes tt c [see panel 
(a) of Fig. 2]. The time of r r (t) = T t (t), where 
adiabatic-diabatic transition occurs, corresponds to the 
freeze-out time in KZ theory JlB|. In the single-stable re- 
gion, tt > tt c , introducing tt(t) — tt(t) — tt c , we have 

2f2(A+Sl c ) 

— n a (2h+n c ) Tq ' 



and t 4 
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where t = l/tt c . For slow transitions, r g ^> 1, we find 
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In the bi-stable region, f2 < f2 c , introducing = 
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For slow transitions, it is easy to findji 
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z = 1 and = 1/2, for continuous quantum phase tran- 
sitions @,[ioj|. 

Due to the adiabatic condition becoming invalid near 
the critical point, the defect modes could be stimulated. 
In Fig. 2, we show the dynamics of SBTs in the MF 
system. In our simulation, the initial state is chosen as 
a symmetry-broken ground state in the bi-stable region, 
and tt is ramped up from below to above tt c . For different 



ramping rates, we calculate the longitudinal polarization 
and the fidelity to instantaneous ground states. The re- 
sults show that the MF defect modes appear as damped 
oscillations with /^-dependent amplitudes [see panels (b) 
and (c) of Fig. 2]. The defect generation via dynamical 
insatiability is similar to the vortex nucleation in rotating 
BECs [2lj |. Interestingly, the defect modes become sig- 
nificant after the revival of adiabaticity when the system 
has passed the critical point, but not after the breakdown 
of adiabaticity when the system approaches the critical 
point. This indicates that the system has no sufficient 
time to complete a whole intrinsic oscillation in the in- 
terval between the breakdown and revival of adiabaticity. 
Because of this, the system tries to keep its present state 
when the adiabatic condition is invalid. 
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Figure 2: The MF dynamics of SBTs in Hamiltonian (2). 
(a) The Tt and r r , (b) the longitudinal polarization (Nf — 
Ni)/2, and (c) the fidelity to instantaneous ground states 
versus ft/Ec for different values of (3. Parameters are chosen 
as N = 100, G u = G TT = 1.0, G n = 2.0, and E (n = E or . So 
that, E c = -2.0 and A = 0. 

To explore the correspondence between MF and FQ 
dynamics of SBTs, we compare the dynamics of Hamil- 
tonians (1) and (2). In Fig. 3, we show the FQ dynam- 
ics corresponding to Fig. 2. In the FQ system, simi- 
lar defect modes appear as damped oscillations. How- 
ever, their critical points and oscillation regimes are dif- 
ferent from the MF ones. First, the FQ defect modes 
appear after the critical point between quasi-degeneracy 
and non-degeneracy, which is not identical to the MF 
critical point between bi-stability and single-stability [see 
panel (b) of Fig. 3]. Additionally, the FQ defect modes 
have /?- independent amplitudes [see panel (c) of Fig. 3]. 
Furthermore, starting from a ground state in the single- 
stable region, the MF defect modes will always appear if 
the system passes its critical point. However, if is suf- 
ficiently small, there is no FQ defect mode. That is, the 
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FQ system always remains in its instantaneous ground 
state. This means that, the MF dynamics well coincides 
with the FQ one before the system reaches the first crit- 
ical point and then it breaks down. The MF breakdown 
occurs at its FQ critical point if the system goes from 
bi-stability to single-stability. In contrat, the MF break- 
down occurs at its MF critical point if the system goes 
from single-stability to bi-stability. This anomalous MF 
breakdown dependent on approaching direction differs 
from the conventional MF breakdown related to inter- 
action strength and total particle number. 



down induced by SBTs depends on the approaching di- 
rection. In a transition from the normal to polarized 
regions, the MF breakdown occurs at the MF critical 
point. However, in a transition from the polarized to 
normal regions, the MF breakdown occurs at the FQ crit- 
ical point. This anomalous MF breakdown broadens our 
conventional understanding for MF breakdown related to 
interaction strength and total particle number. 
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Figure 3: The energy spectrum and FQ dynamics of SBTs 
shown in Fig. 2. (a) Energy spectrum of four lowest states 
of Hamiltonian (1), (b) and (c) quantum expectations of the 
longitudinal polarization (Nf — Ni)/2. 

To observe the dynamical SBTs, one should prepare a 
two-component BEC jjjj]. The negative Ec with all posi- 
tive s-wave scattering lengths can be obtained via Fesh- 
bach resonances. To avoid the phase separation induced 
by strong inter-component repulsion, one has to trap the 
BEC within a sufficiently strong potential. The inter- 
component coupling can be realized by Raman and/or 
radio-frequency fields, and then the coupling strength is 
controlled by the field intensity. By ramping the field 
intensity up or down, one can push the system pass its 
critical points. To detect the longitudinal polarization, 
one should count the atoms in each component via state- 
dependent fluorescence or spatial imaging. 

In conclusion, we have studied both MF and FQ dy- 
namics of SBTs in a coupled two-component BEC. We 
analytically obtain universal KZ scalings for MF dynam- 
ics and numerically explore the correspondence between 
MF and FQ dynamics. The MF dynamics well coin- 
cides with the FQ one before reaching the first critical 
point and then it breaks down. Because the FQ criti- 
cal point is not identical to the MF one, the MF break- 
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